Currents and pseudomagnetic fields in strained graphene rings 
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We study the effects of strain on the electronic properties and persistent current characteristics 
of a graphene ring using the Dirac representation. For a slightly deformed graphene ring flake, one 
obtains sizable pseudomagnetic (gauge) fields that may effectively reduce or enhance locally the 
applied magnetic flux through the ring. Flux-induced persistent currents in a flat ring have full 
rotational symmetry throughout the structure; in contrast, we show that currents in the presence of 
a circularly symmetric deformation are strongly inhomogeneous, due to the underlying symmetries 
of graphene. This result illustrates the inherent competition between the 'real' magnetic field and 
the 'pseudo' field arising from strains, and suggest an alternative way to probe the strength and 
symmetries of pseudomagnetic fields on graphene systems. 

PACS numbers: 73.22.Pr, 73.23.Ra, 61.48. Gh, 73.23.-b 



The appearance of gauge fields in graphene is a beau- 
tiful and experimentally accessible example of a situ- 
ation where concepts of condensed matter and quan- 
tum field theory converge on a physical systemP^ Ex- 
perimental evidence of 'bubble' formation on particu- 
lar graphene growth processes^ and contr ollab le routes 
to manipulate graphene bubble morphology^' have mo- 
tivated numerous works addressing different aspects of 
these effects. In particular, the theoretical description 
of strained graphene has been developed significantly, 
exploring how its electronic properties are affected on 
curved and strained surfaces. It is known that elastically 
deformed graphene can be mapped onto the Dirac for- 
malism in the continuum limit by including pseudomag- 
netic fields^H^ and Fermi velocity renormalizationP giv- 
ing rise to local quantities that depend on strain but do 
not break time reversal symmetry. A recent contribution 
reported a space-dependent Fermi velocity leading to in- 
teresting experimental consequences.^ Lattice-corrected 
strain induced vector potentials in graphene have a lso 
been discussed within a tight-binding scenario pJil2l a i_ 
though these corrections do not c ontribute to the pseu- 
domagnetic field distribution! 12 ! 13 ! 

Similar to other confined systems with periodic 
boundary conditions, magnetic flux-dependent persis- 
tent currents and conductance oscillations are ex- 
pected for graphene rings in an Aharonov-Bohm (AB) 
geometry! 14 " 16 ! Several experiments have verified the 
presence of AB conductance oscillations with different 
visibility for different device geometries G3HH1 A recent re- 
view of quantum interference in graphene rings discusses 
open questions in the field! 2 ^! Interestingly, the 'infinite 
mass' confinemenlPSt2ll that requires null current density 
across the boundaries results in persistent currents that 
are 'valley polarized' in the presence of magnetic fhix™=! 
suggesting that graphene quantum rings would be an ex- 
cellent system to analyze the effects of induced curva- 



ture. The main result of the present paper is indeed to 
show that while a flat (unstrained) graphene ring in the 
AB geometry sustains persistent currents with full rota- 
tional symmetry, unavoidable strains in typical systems 
would result in inhomogeneous distributions of currents. 
In other words, while the strains alone would result in 
zero net persistent current (since time reversal symmetry 
is not broken by the pseudo magnetic fields) , the compe- 
tition with the AB flux induces spatially inhomogeneous 
current distributions on the system. This effect can be 
seen to arise both from a local rescaling of the Fermi ve- 
locity as well as by the appearance of gauge fields that 
result from the elastic deformations. As such, the persis- 
tent currents originated by the magnetic AB flux acquire 
a local character that follows the strain fields. Moreover, 
as the corresponding length scales of the current inhomo- 
geneities are given by the strain fields, one can imagine 
using this effect to measure the strain distribution via 
a scanning magnetometer that would be sensitive to the 
induced currents. Alternatively, properly designed strain 
fields would be used to produce desired current patterns. 

Strained graphene. Within a tight-binding model, the 
effects of lattice deformations may be incorporated into 
the hopping integrals t n between nearest-neighbors p3 
so that t n = t + 8t n = t (l — PeijS^S^/a 2 ), where 

(3 = |<91ogio/(9loga| « 3 in graphene, S n are the nearest- 
neighbor vectors of a given atom at lattice site n, and 
to and a are the nearest-neighbor hopping integral and 
distance in the unstrained system, respectively. Indices 
i and j represent directions on the 2D plane, with re- 
peated index summation convention throughout. The 
strain tensor, = | (djUi + diUj + dihdjh), is charac- 
terized by Uj and h, the in- and out-of-plane deforma- 
tions, respectively.^ 

The resulting Hamiltonian in the presence of inho- 
mogeneous strain and external magnetic field given by 
B = V x A ext , can be written in a generalized Dirac 



2 



form (in the K valley) as 

H = -iv kj a k (dj + i^Af ') -iv o a T j+v^a^f , (I) 

where <7j are Pauli matrices, and the vector potential 
with trigonal symmetry arising from strains is given by^ 

%\V = p-(^r)(e**-evv-**y) > (2) 
Z7T y la J 

with $o = h/e. The renormalized Fermi velocity tensor 
in Eq. (TTT) is position-dependent 



Vij(r) = v [I, 



(2ey + r\ijd kk ) 



(3) 



where I is the 2x2 identity matrix and vq — 3ato/2 w 
10 6 m/s (with h = 1). Inhomogeneous strains also yield 
an effective geometric vector potential 



2 Utt n 



(4) 



At the K' valley, both and Ti(r) are the same, 

while the vector potential A s Mr) changes sign, preserv- 
ing overall time-reversal symmetry of the system in the 
absence of A ext . 

Diagonalization of the Hamiltonian in a ring geometry 
results in interesting eigenstates and persistent current 
patterns for states at and near the Dirac point (charge 
neutrality point). Wc first summarize the results for a flat 
graphene ring in a magnetic fluxp' to provide a suitable 
framework for the ring with deformation. 

Unstrained graphene ring. We consider a 'flat' ring 
threaded by a magnetic flux <£>, with A ext — ($>/2nr)6, 
so that the Hamiltonian is given by 



Ho 



-iv 



Ai(0)3 r + A 2 (0) 







(d e - 







(5) 



with Ai (9) = <J X cos 6 + <7 y sin 9 and A 2 (9) = —<J X sin ( 
a y cos 9. The wave functions for energy E are^ 



ise ie 4>rn+i{r) 



(6) 



with 4> rTl (r) = AfnJm(kr) + BmYm(kr), where m = 
rn + m = 0, ±1, ±2, ... is the orbital angular mo- 

mentum, Jfh and Ym are Bessel functions of first and 
second kind, respectively, k = \E\/Jwq, and s = sgn(-E). 
The upper (lower) component of the spinor corresponds 
to ipA ("0s) m Eq. ^ for the K valley. 

The energy spectrum is obtained from the transcen- 
dental equation that arise s after imposing infinite-mass 
boundary conditions JliJSiJ gj ven by ip B (r J 9)/ipA(r,9) — 
irn ■ re lT ° , where the normal to the boundaries is n = ±r 
for the inner (— ) and outer radius (+) of the ring. The 
K and K' valley eigenstates for a given m have also 
a total angular momentum j, defined for the operator 



J z = —idef, + <t 2 t/2, where r = ±1 identifies the valleys; 
K and K' states with opposite momentum j' = —j are 
related by m! — — (m + 1). Notice that the boundary 
conditions do not mix the valleys and i n fact break the 
valley degeneracy for non-zero flux! 14 ! 24 ] 

Figure [lja) shows the energy spectrum of a graphene 
ring vs. magnetic flux; dashed and continuous (blue) lines 
(upper part of the graph) indicate results for K and K 1 
valleys with m and m! values, respectively. Notice the 
quadratic dependence on $ of these levels, which breaks 
valley degeneracies in general. The figure also shows the 
lowest energy levels for a deformed ring (thicker red lines, 
lower part of graph), to be discussed later. 
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FIG. 1. (Color online) Unstrained graphene ring, (a) Energy 
spectra vs. magnetic flux for ring with internal and external 
radii given by r\ = 50 and r^ = lOOnm, for different quantum 
states: m (m' = — (m + 1)) integer denotes results for K 
(K') valley given by dashed (continuous) lines. [Thicker (red) 
lines near E ~ 17meV show lower levels for the same ring, 
but with a Gaussian deformation.] (b) Flat ring eigenstate 
with m = (K valley); electronic probability distribution 
along the ring for $ = 0, and (c) $/$o = 6. Main panels 
show the two spinor components separately; insets show total 
distribution, \iI>a\ 2 + \iPb\ 2 - On right column, corresponding 
current densities for K valley, (d) Notice counter-propagating 
edge currents for $ = 0, indicated by red arrows, evolve to a 
current distribution mostly on the outer radius for large $ in 
(e). 

The charge and current densities that satisfy the con- 
tinuity equation for unstrained graphene are given by 
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p = and Jj — (vo/hjt/j^Ta-jil;^ Typical results for 
the spinor components \^a\ 2 and l^sl 2 for m = are 
depicted in Fig. [T] along the radial direction, for both 
(b) $/$o = 0, and (c) 6. As expected, increasing flux 
causes the charge density to be driven to the outer edge 
of the ring as the energy of the state increases. (Notice 
the m = state is not in general the lowest state as <I> in- 
creases.) Also shown in Fig.JTJd) and (e) are the current 
probability densities, highlighting the strong dependence 
on the magnetic flux. Notice that the current for valley 
K at zero flux is given by nearly compensated counter 
propagating edges, while as the flux increases, the cur- 
rent in the inner edge disappears. As we will see below, 
the interaction with the pseudomagnetic field generated 
by the deformation gives rise to a much intricate current 
pattern. 

Strained graphene rings. We now consider an out of 
plane deformation given by a circularly symmetric Gaus- 
sian shaped described by h — Ae~ r l b . The strain ten- 
sor is then 



af(r) 



sin v cos ( 



sm o cos ( 
sin 2 6» 



(7) 



where f(r) = 2 (r 2 /b 2 ) e- 2r2 / b \ with a = A 2 /b 2 char- 
acterizing the strength of the strain perturbation. The 
strain is inhomogeneous and, as a consequence, the 
geometric gauge field T is nonzero and the renormal- 
ized Fermi velocity changes along the ring. The space- 
dependent velocity is given by, 



f'o 



-/3a 



fir) 



I+-R{29)a z 



(8) 



where R (9) = I cos 9 — ia y sin# is the rotation matrix 
through an angle 9 in the counterclockwise direction. The 
resulting Dirac cone becomes elliptical, with radial and 
angular components 



Wo 



-3(3a 



f(r) and v e = v 



-/3a 



while the gauge fields are 



(9) 



r. r = 



-/3a 



3r 



and Ta = 



(10) 



These changes can be seen as perturbations of the Hamil- 
tonian in Eq. ([5| given by 



-iv 



-/3a 



f(r) 



A 1 (9)d r + A 2 (9)- 



(11) 



with d r = §<9 r + I - §5 and dg = \ [d e + i^j . The vec- 
tor potential perturbation V<i — vq (y-^f^ /(r)Ai(— 28) 



is associated with a pseudomagnetic field Bgt, given by 

(12) 



The eigenvalue problem with H = H q + V\ + Vi can be 
solved using perturbation theory on the parameter a up 
to second order, keeping sufficient states to achieve full 
convergence of the results. 

We now present our main results for strained rings con- 
sidering the Gaussian perturbation with characteristic 
system parameters: A — 7nm and b = 70nm, with a rela- 
tive deformation a = 1%, and the ring radii used in Fig. 
[lja). The two lowest states of the spectrum corrected 
by the Gaussian deformation are shown in Fig. [lja) in 
thick (red) curves near the bottom of the panel, both for 
solutions near K (dashed lines) and K' (solid) valleys. 
We find that the main correction comes from the V2 per- 
turbation which contains the effects of the strain-induced 
pseudo magnetic field, and produces energy shifts for the 
ground state as high as 10%. 




FIG. 2. (Color online) Deformed graphene ring. Contour 
plot of the local density \4>a\ 2 + \ tps\ 2 for m — state at (a) 
$/$o = 0, and (b) $ = —0.5. A negative flux, as shown, 
pushes density towards the inner ring, while a positive flux 
would shift weights to the outer radius, (c) Current distri- 
bution for $ = and m — 0; for the K valley, (d) Current 

variation, A J = J — Jflat, (if -valley) where Jflat is the cur- 
rent for a flat ring-shown in Fig. [TJd)-can be seen as the net 
effect produced by the deformation. Notice six vortices with 
alternating circulation and cores near maxima/minima of the 
pseudomagnetic field (see Fig. [3ja)). 

The pseudomagnetic field produced by the Gaussian 
deformation in this system (see Fig. [3]ja)) has the under- 
lying trigonal symmetry of the graphene lattice!^ In this 
case, the field amplitudes reach ~ 1.2T, and, as expected, 
when averaged over the entire ring, the net pseudo field 
vanishes. 

When both the external magnetic flux and deformation 
strains are considered, the superposition of fields with 
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different symmetry strongly affects the electronic states 
and induces inhomogeneities in the probability density 
distribution. Typical changes in the spatial pattern of 
electron density are shown in Fig. [2ja) and (b), where 
zero and finite fluxes are considered. The local density 
is shown by colored projections along the ring. A finite 
flux, either positive or negative, enhances the amplitude 
modulations of the local probability density seen for zero 
flux: a negative flux shifts the maxima towards the inner 
radius of the ring, while a positive flux pushes the den- 
sity towards the outer radius, as one would expect from 
classical Lorentz force considerations. 

The current density J over the strained graphene ring 
is displayed in Fig. [2jc) , for the lowest state in the K val- 
ley. The current density trends are represented by a set 
of small (red) arrows, revealing intricate current configu- 
rations. The current density exhibits local maxima with 
trigonal symmetry, centered near regions of largest varia- 
tion in the local density. The last panel, Fig.[2jd), shows 
the probability current variation, A J = J — Jflatj where 
Jfiat is the persistent current in the flat or unstrained 
ring. Notice that A J exhibits six vortices with alternat- 
ing circulation and cores centered on regions of extremal 
values (positive or negative) of the pseudomagnetic field 
in Eq. @, Fig. [fa). 

The spatial modulation of the persistent current in the 
ring can also be seen from the angular profiles shown in 
Fig. |3jb) for different flux values. The current density 
achieves maximal values along nodal lines of the pseudo- 
magnetic field distribution. Similarly, local minima are 
related to regions of maximum pseudo field amplitude, 
both in positive and negative directions-see also polar 
plots in Fig. [3]ja). As the pseudo field does not break 
time-reversal symmetry, the net current must be zero in 
the absence of external magnetic flux. This indeed hap- 
pens when considering the contribution of the other val- 
ley (K' , coming from the state with m' = — 1), which has 
fully inverted symmetry in 9. 

We analyze the role of the pseudomagnetic field con- 
tribution to the total current (taking both valleys into 
account) by looking at different angular values along the 
ring. Figure^c) shows the total current at 9 — n/6, 7r/3, 
and 7r/2, in comparison with the (angle integrated) cur- 
rent for both the unstrained and deformed graphene ring. 
All curves present the expected sawtooth behavior with 
flux. However, the slope of the curve and the value near 
zero flux are clearly angle and strain dependent. Notice 
in particular the jump reversal near $ = for 9 = it/ 6 
and 7r/3, associated with the circulation around the vor- 
tex at 9 ~ 7r/6. The competition between the external 
magnetic field and the pseudo field not only results in 
inhomogeneous current distributions with vortices, but 
also in very different total current dependence with flux 
$. Notice that the strain decreases the total current dis- 
continuity near $ = 0. 

Conclusions. We have shown that the strain effects 
arising from a Gaussian 'bubble' deformation of the 
graphene ring result in a distribution of pseudomagnetic 
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FIG. 3. (Color online) (a) Map of pseudomagnetic field in 
Eq. ( |12[ ). Traces show polar plots of the current with $/$o = 
—0.5 for ring without (circular) and with (star) strain, (b) 
Angular profiles of current density through the ring for the 
lowest K valley state and different fluxes (A$/$o = 0.1), 
from $ = to $/$o = —0.5, as indicated by the arrow. 
Dashed flat line near bottom shows current density for un- 
strained ring at $ = 0. Notice mean value of current in- 
creases with strain and flux, (c) Flux dependence of the cur- 
rent through the ring for the lowest state: m = 0, -valley 
from "l?/$o = —0.5 to 0, and m = —1 for Jf'-valley from 
$/$o = to 0.5, at different angles along the ring. Curves 
for strained rings are shifted up for clarity; dashed lines in- 
dicate zero current in each case. Bottom traces show angle- 
integrated current for both strained (solid) and unstrained 
(dotted) rings; notice similar slopes, although much smaller 
discontinuity near $ = with strain. 



(gauge) fields that have trigonal symmetry, in agreement 
with the underlying symmetries of graphene. While the 
currents induced by these pseudo-fields would identically 
vanish, an external magnetic flux makes possible the ob- 
servation of the full spatial distribution of currents due to 
strain. As a result, strain fields change the nature of the 
ground state and modify the amount of current present 
in the device. 

Our discussion has focused on the infinite-mass bound- 
ary condition. We find also that the zig-zag boundary 
condition, which does not mix valleys either and shows 
a completely different spectrum, yields qualitatively sim- 
ilar results to those presented here. It is clear that an 
experiment would typically exhibit a more complex edge. 
However, because the effects are produced by strains 
fields, there will always be inhomogeneous current dis- 
tributions upon the application of a flux. 

Furthermore, other geometric structures with strong 
strain fields would also produce complex patterns of in- 
duced currents reflecting the pseudomagnetic field dis- 
tribution. These could be explored locally by applying 



5 



a weak magnetic field, to play the role of the external 
flux, assuming that the scanning current measurement 
device has spatial resolution better than the characteris- 
tic length scales of the gauge field distribution. A scan- 
ning magnetometer appears as one of the ideal instru- 
ments to reveal these effects. 
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